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Abstract 


Four configurations of MLS Flare Elevation Systems (E^) that can be 
considered reasonable and practical in actual implementation are identified. 
Each of these are analyzed and compared with respect to (a) computational 
requirement, (b) required coverage, and (c) accuracy including altitude and 
sink-rate estimation error performance. 


I. Configurations 

The four (4) configurations are based on the coordinates (conical or 
planar) and the orientation of the beam center (conventional or nonconven- 
tional; parallel or perpendicular to the runway center line, respectively). 

(i) Conical/conventional (the most likely U.S, candidate) 

(ii) Conical/nonconventional 
(iil) Planar/ conventional 

(iv) Planar /nonconventional (ALI; Australian candidate). 


II . . Computational (or Algorithm) Requirements 

By the computational requirement is meant the difficulty or the sim- 
plicity of obtaining the altitude from El^, DME and Azimuth information# 

The following notations are used. 


r 

¥ 

e 

x,y,z 


= DME slant range 
= Azimuth Angle 
= EI 2 flare' elevation angle 

“ Runway referenced rectangular coordinates of the 


aircraft 

XAZ, ZAZ =» Az/DME (colocated) antenna site. 

XE^, YE1 2> z E 1 2 " El 2 antenna site 


(A) MLS Equations from x,y and z 

The range (DME) and the azimuth (conical) equations are the same 
regardless of the configurations: 

DME Equation: 

(1) r » [ (x-XAZ) 2 + y 2 + (z-ZAZ) 2 ] 1 / 2 

(2) Az Equation: 

* = sin" 1 (— ) 

The El 9 equations are not the same for the various configuration 
They are given next for each configuration. 
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(1) Conical/conventional 


(3) 


-1 r z-ZEl? 

G “ ^(x-XE^)^ + (y-YEl^* + (z-ZEl,,)*) 


(ii) 

Conlcal/nonconventlonal 



Same as (3) 


(Hi) 

Planar /conventional 

(4) 


E tan 1 x -XSl 2 * 


(Iv) 

Planar/nonconventional * 


(5) « - tan- 1 ( } 


NOTE: x, y, z are expressed in the runway center line coordinate system 

shown below: 



(B) Altitude Equation from MLS 
(i) Conical/conventional 


For this configuration Equations (1, 2, and 3) must be 
used to determine the altitude precisely. The solution involves 
a quadratic equation: 

(6> Z - ZE1 2 - i(K-XEl 2 )Z + (y -YEt 2 )2)} 1/2 , 

where 

(7) y = r sin i/> 

and x is a of a quadratic equation 

(8) Ax 2 + Bx + C = 0 


with 
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A 


1 

1- tan' 1 


El 


(9) 


B = -2 


{ XAZ + XE1 2 


tan 2 s 
l-tan z eT 

.2 - vt?i2 


•> 


9 9 . . , tan 2 e r XEl 2 + [r SIN t \> -YEl?] 2 ^ 

-r z cos z i{/ + XAZ + f : t an 2 ~7 { 2 > 


Equations (6) through (9) are exact solutions in the sense 
that there is no approximation involved except the antenna heigf 
was assumed the same. 

For the same order of exactness the following recursive 
(or iterative) algorithm can be used: 

+ 

Assume an estimate of altitude, Z . 

Obtain v = r sin iji 

Using Z+ & y, obtain x from Equation (1). 

Using z" 1 " & y & x and e, obtain z from Equation (6) 
Use z as the estimate and go back to step (ii) 
and repeat. 

NOTE: This recursive algorithm has been evaluated for the compu- 

tation of altitude from DME, Az and El^ on the glide-slope, 
and it converged rapidly. 

Algorithm (10) still involves considerable computation, 
therefore the planar approximation 

(II) Z = ZE1 2 - (r - XDIST) tan e 

is used, where 

XDIST = distance between DME and El 2> 

As will be shown later, this approximation turns out to be 
surprisingly good in spite of its simplicity. 


Step. (i) 
Step (ii) 

(10) Step (iii) 

Step (iv) 
Step (v) 


(ii) Conical/nonconventional 


In order to determine altitude with this configuration, 
either Equations (6) through (9) or the recursive algorithm (10) 
must be solved. The planar approximation would not yield good 
results, since the distance to the aircraft is much shorter. 


(iii) Planar/conventional 

The exact solution is given by: 
(12) Z = ZE1 Z + (x - XE1 2 ). tan e 




where 
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(13) 

(14) 

(15) 


(16) 


(17) 


(18) 


y ■ r sin i|» 
and x is a solution of 
Ax 2 + Bx +C “ o, 

« 1 + tan 2 e 

B = -2 [XAZ - XE1 2 • tan 2 e] 

C - XAZ 2 + XEl 2 tan -r 2 cos 2 <p 
2 

The recursive algorithm corresponding to (10) can be derived 
for this case also; however, the computation savings is nil. 

Similar to Equation (11), we can obtain an approximation: 


Z = ZE1 2 - (r - XDIST) . tan e 

This approximation (using DME directly instead of computed 
as in (12)) should yield excellent accuracy. 


(iv) Planar /nonconventlonal 


For this configuration the exact solution is the easiest and 
is given by: 

Z = ZEl 2 -(r ■ Sin i|» - YE1 2 ) • tan e, 

If the aircraft is assumed on the extended runway center- 
line the equation simplifies to: 

Z = ZE1 2 - YE1 2 • tan e, 

since tji a 0. 


Conclusions for Section II 

(1) For conical coordinates (Conf. (i) & (il)) the exact altitude 
determination is very involved. This can be simplified somewhat using 
recursive formulation with very good accuracy; however, the algorithm 
is still involved. For configuration (i) the planar approximation can 
be used with satisfactory results. 


Ill . Required Coverage 

The El 2 coverage requirements are specified by RTCA DO-148 (page 8). The 
Flare signal must be available near the runway surface throughout the touch- 
down zone and shall extend to a range of 0.5 nm from the threshold. For a 3° 
glideslope a nominal flare initiation altitude of 50 feet at 200 feet/second 
landing speed, the above requirement gives 100 feet (approximately 10 seconds) 
of margin to acquire and settle the filter prior to initiating the flare 
maneuver. (From the coverage point of view, the portion from 150 ft. ^ 50 ft. 
altitude is the only advantage the MLS derived flare has over radio altimeter 
since the latter may not be reliable ahead of the threshold due to irregular 
terrain at some airfields.) 



rKUNHW THRESHOLD 



MINIMUM COUP PAGE 
AREA FOR MLS 

el z, function 

APPROACH 

D/EECTtON 


Fig. 1. Minimum Coverage Area for Flare Function 


(A) Configuration (i) 

Assuming the flare antenna is 2500* behind the GPIP toward the 
runway stop-end and offset 450' laterally from the runway centerline, 
the coverage for configuration (i) is given by (see Figure 2) 

lateral coverage = tan ^ , 450' + 100*. = 47.7° 

* 500’ ; 

vertical coverage = tan . 150 ' . = 1.4° 

'“3000' + 3500 y) 


NOTE: The 1.4° vertical coverage is substantially less than 

specified in DO -148 (page IIIA-11), 


(B) Configuration (iv) 

For Configuration (iv) the lateral as well as the vertical 
coverage are more sensitive to the offset distance since (see Figure 3) 


lateral coverage = tan 


-1 


0 


3000' 


D offset - 100 


") 


-1 


vertical coverage = tan ^ 


150' 


D offset - 100 


r) 


where 


D offset = Flare antenna lateral offset distance from the runway 
centerline. 
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The following table summarizes the coverages for different values 
of D offset. 


D offset 
(ft) 

Lateral 

Coverage 

(deg) 

Lateral 

Coverage 

(deg) 

500 

82.4 

20.6 

750 

77.8 

13.0 

1000 

73.3 

9.5 

1250 

69.1 

7.5 


Table 1. Lateral & Vertical Coverage for 

Configuration (iv) 


IV . Accuracy and Filtering Performance of Configuration (i) and (iv) 

In this section the following cases are examined for a nominal expo- 
nential flare at 200 ft/sec landing speed and 3° glide slope. (See Figure 4.) 

(1) Configuration (1) 

(l.a) Planar approximation (conicity) altitude error 

(l.b) Altitude and sink-rate estimation errors with complementary 
filter with planar approximation. 

(l.c) Same as (l.b) except with noncomplementary filter (i.e., no 
acceleration information.) 

(2) Configuration (iv) 

(2. a) Altitude and sink-rate estimation error with complementary 
filter and zero lateral (y) deviations 

(2,b) Same as (2. a) except with noncomplementary filter. 

Computer runs were made with the conditions and equations shown below. 
(Only one run was made for each set of conditions; however, the same noise 
sequence was used for each configuration so that comparitive results were 
obtained.) 

Constant lateral (y) deviations were employed to approximate the effects 
of aircraft lateral errors on measurement performance and to simulate the 
coupling between the aircraft lateral and vertical motions. The magniture of 
the aircraft high frequency lateral motion is generally small and hence the 
induced high frequency vertical errors will also be small; however, the low 
frequency lateral motion induces a slowly varying vertical path following error 
which cannot be filtered effectively. This coupling effect will be most pro- 
nounced with configuration (iv) with the simplified altitude computation; 
however, the problem is substantially reduced if the exact altitude computation 
is utilized. 





S'c/C r 7£' ^<7/7 £> €>3G 2T O/VX 00^-79^7 (g) 

/'SJ.te/ ^LH9/7=r 7iX/JL/VSf\/Q=/X9 7t/A//WOA/ > S9//)7>/-/ 
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Following conditions were ut-ed. 

°E1 2 Ant. site for Conf. (i) 

XE1 2 - 2500', YE1 2 - 400', ZE1 2 » -10*. 
°E1 2 Ant. Site for Conf, (iv) 

XE1 2 « 0, YE1 2 = 500', zei 2 = -10‘. 

°E1 2 Errors ** 

bias = 0.015, 0., - 0.015° (+ lo) 

•Noise - 0.016° (la) 

°DME Error 
bias = 0 

Noise « 20 St,, (la) 

Equations 


0 Nominal Exponential Flare Flight-path 


x . , « x + x A 
n+1 n 


^n+1 = y n = const a ”30 ,, 0, 30' 

foi 

flight path) 


Z n+ ^ = Z^ + Z^ A (Z n = const for 3° glide slope portion of 


Z jn+1 " TtZ T.D. “ Z n ] 
Z n+1 = exp (-A/t) • Zn 

Z n+1 = " Z n+1 /T 


T - Z INT /(Z TD - Z INT ) 


Exponential Flare 
Flight-path portion 


x D = -1431 ft, x = 200 ft/sec, Z 0 = -75 ft, Z 0 = 10.5 ft/sec, 

Z INT “ “ 50 ft ’ ^XNT = 10,5 ft/,sec » ^TD = 2,3 ft / sec > 
t =* 6.05 sec. 


**Note that since the El-2 antenna location for configuration (iv) is closer 
to the touchdown point than for configuration (i), the configuration (iv) angular 
error can be two to three times larger for the same linear altitude error. 
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0 Altitude computation for Conf. (i); see Equation (11) § II 

Z' a ZE1, - (r - 11355.0) tan e 
n 2 ' n n 

0 Altitude computation for conf. (lv), see Equation (18), § V 

Z* a ZE1_ - YE1„ • tan e 
n L l n 

where r is the DME measurement and e is the El„ measurement. Roth 
n n / 

contain the errors given previously. 


Filter Equations 


Complementary and noncomplementary filters were implemented in 
sequential versions, since the sequential format is much simpler with 
a digital computer. 


(a) Complementary Filter 



Acceleration 

Complementation 

Term 



(b) Noncomplementary Filter 



P 

( ) n - a-priori estimate, ( ) = a-posterrori estimate, 

( )’ = observed or measured value at time n. 
n 

g cl “ 0, ° 81 ’ S c2 = 0-017, gn^ “ 0.156, gn 2 => 


0.066 



Results 


The results of computer runs are summarized below. The heading 
numbers correspond to the cases itemized at the beginning of this 
section. 

(1) Configuration (i) ; Conical/C^nventlonal (See Tables 2a and V, 

(l.a) The planar (Equation 11) approximation was found to be 
surprisingly accurate. Along the nominal exponential flight-path and on 
the azimuth" £. the computed altitude error was between -0.34 feet at 
75' altitude to +0.34 feet at T.D. Since this configuration (and the 
approximation) is insensitive to lateral y-deviations, the result will 
hold for reasonable y-deviations, say within + 30 feet. 

(l.b) The maximum altitude estimation error was less than 1.6 feet 
(Iff) and the sink-rate estimation error was less than 0.3 ft/sec (lo). 

The estimates were very smooth. The bias in El, was the major source 
of error. Even though not shown, DME noise of 20 ft (10) contributed 
very little error. The apparent touch down distance (where the filter 
indicates the aircraft touched down as opposed to where it actually 
touched down were nicely clustered. No noticeable y-deviation effects 
were observed. 

r 1 . c ) The performance for the noncomplementary filter was decidedly 
inferior compared to the complementary filter case. The maximum 
altitude error was s; 2.5 ft. (10) and the maximum sink-rate error ~ 

2.6 ft/sec. The estimates were very noisy. This is to be expected 
since the filter's band-width is eight (8) times that of the comple- 
mentary filter, and thus allows more MLS noise to pass through. (How- 
ever, the band-width may not be reduced arbitrarily, since the dynamic 
delay error will begin to dominate.) The apparent touch down distances 
were more dispersed and the bias in EI 2 measurements had more prominent 
affect on the errors. 

(2) Configuration (iv) ; Planar/nonconventional (Set Tables 3a and b) 

(2. a) The maximum altitude estimation error was less than 3.3 ft. and 
the sink-rate error less than 0.5 ft/sec. with complementary filtering. 

The estimates were very smooth. The apparent touch down points were 
scattered compared to (l.b). The bias and the noise of EI 2 had much 
less affect than in (l.b), since linear errors due to El, error are 
much smaller due to shorter distance to the aircraft. The major source 
of error was due to the neglected lateral y-deviation. 

(2.b) Without complementary filtering the maximum altitude estimation 
error was less than 2.7 ft. and the sink-rate error less than 1.63 ft/sec. 
The estimates were somewhat noisier than the (2. a) case, but smoother 
than the (l.c) case. This is because the errors in El 2 do not contribute 
greatly to the altitude computation. The apparent touch down points were 
similar to the (2. a) case. The largest errors occurred at or immediately 
after the flare initiation, since this is where the vertical acceleration 
is largest and hence the most severe maneuver dynamics effect takes 
place. Soon after the filter's performance approaches that of (2. a). 

The affect of y deviation on vertical error was essentially the same 
as shown in (2. a). 
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y- dev. ■ 
* (feet) 

E12 bias 
(deg) 

(o; 

Alt. Error 
@ 50’ 

(ft) 

CD 

Max Alt. 
Est. Error 
(ft) 

(2) 

Max Sink 
rate 
error 
(ft/sec) 

(3) 

Apparent 

T.D."' 

(ft. from 
GPIP) 

1" 

30 

0.0 

0.8 

0.5 

0.08 

946 


0 *' 

0.00 

0.8 

0.5* 

0.04 

946 


-30 

0.0 

0.7 

0.6 

0.09 

946 


30 

0.015 

1.7 

1.1 

0.16 

986 


0 

0.015 

1.7 

1.0 

0.16 

986 


-30 

0.015 

1.7 

0.9 

0.15 

976 


30 

-0.015 

-1.1 

1.5 

0.21 

906 

0 

* 

-0.015 

0.1 

1.5 

0.22 

906 

-30 

-0.015 

-0.2 

1.6 

0.22 

906 


(0) Alt Error at 50' = Alt act Alt comp at 50' Alt 

(1) Maximum estimated altitude error - Alt act - Alt est during flare 

(2) Maximum estimated sink-rate error = (Alt) est during flare 

(3) Apparent' longitudinal T.D. point - where the filter thinks the a/c 
T,D.ed» i.e. , '.'here tha Alt est = 0.0 

H 

Table 2. a) Conical/Conventional with 

Complementary filter (t =* 2.36 sec) & 

Planar Approximation 

^Nominal value for T.D. - 946’ 


ORIGINAL PAGE IS 
OP POOR QUALITY 
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y-dev. ‘ 
(feet) 

EX2 bias 
; (deg) 

(0) 

Alt Error 
@ 50' 
(ft) 

(1) 

Max. Alt. 
Est , 

Error ("ft) 

(2) 

Max. Sink 
Rate Est. 
Error (ft/sec) 

(3) 

Apparent 
T.D.* (ft 
from GPIP; 

30 

t 

0.0 

0.8 

1.5 

2.45 

936 

0 

0.0 

0.8 


2.52 

936 

-30 

0.0 

0.7 

1.6 

2.51 

946 

30 

*• 

0.015 

1.7 

2.1 

2.56 

1006 

« 

0 

0.015 

1.7 

2.1 

2.57 

1006 

-30 

0.015 

1.7 

2.0 

2.57 

1016 

30 

-0.015 

-i.i 

2.4 

2.47 

926 

0 

-0.015 

0.1 

2.4 

2.46 

926 

-30 

-0.015 

-0.2 

2.4 

2.46 

926 


(0) Alt. Error at 50' Alt. = Alt act - Alt com p at 50 ' Alt. 

(1) Maximum estimated altitude error = Alt act - Alt est during flare 

(2) Maximum estimated sink-rate error = (\lt) act “ (Alt) est during flare 

(3) Apparent longitudinal T.D. point = where the filter thinks the a/c 
T.D. ed, i.e., where the Alt est =0.0 

Table 2.b) Conical/Conventional w.ith 

Non-Complementary filter (t = 0.295 sec) & 
Planar Approximation 

♦Nominal value for T.D. 3 946' 


































y- dev. 
feet 


E12 bias 
(deg) 



0.0 


0.0 


0.015 


0.015 


0.015 


0.015 


0.015 


•0.015 


(0) 

Alt Error 
0 50' 

(ft) 

(1) 

Max. Alt. 
Est. 

Error (ft) 

-2, 

.3 

2.8 




( 2 ) 

Max. Sink 
Rate Est. 


2 


0.03 


0.48 


0.38 


0.03 


.5 


0.42 


0.04 


0.45 


1006 


946 


886 


1026 


946 


886 


986 


926 


846 


(0) Alt Error at 50' Alt = Alt act - Alt com p at 50* Alt. 

(1) Maximum estimated 'altitude error = Alt act - Alt esC during flare 

• * 

(2) Maximum estimated sink-rate error = (Alt) act - (Alt) est during flare 

(3) Apparent longitudinal T.D. point = where the filter thinks the a/c 

T.D. «d, i.e., where the Alt est =0.0 • • 


Tatie 3. a) Planar/Non-conventional With 

Complementary filter (t - 2.3 sec) 

*Hominal value for T. D. = 946' 


ORIGINAL PAGE IS 
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y-dev. 

(feet) 

E12 bias 
(deg) 

(0) 

Alt Error 
@ 50' 

(ft) 

• (1) 
Max. Alt. 
Est. Error 
(ft) 

(2) 

Max Sink 
Rate Est. 
Error (ft/sec 

(3) 

Apparent 
T.D. "(ft 
) frin GPIP) 

3° 

0.0 

-2.3 

2.3 

0.57 

1006 

0 

0.0 

0.0 

0.3 

1.01 

926 

-30 

0.0 

2.5 

2.5 

1.63 ; 

886 

4 30 

0.015 

-2.2 

2.2 

0.57 

1026 

0 

0.015 

0.2 

0.2 

1.01 

926 

-30 

0.015 

2.8 

2.7 

1.63 

906 

30 

-0.015 

-2.4 

2,5 

0.57 

986 

0 

-0.015 

-0.1 

0.4 

2.0 

926 

-30 

-0.015 

2.3 

2.4 

• 1.63 

886 


(0) Alt Error at 50 1 Alt = Alt acc - Alt corap at 50' Alt 

. ♦» 

(1) Maximum estimated altitude error *» Alt act - Alt during flare 

♦ • 

(2) Maximum estimated sink-rate error = (Alt) acC - (Alt) est during flare 

(3) Apparent longitudinal T.D. point = where the filter thinks the a/c 
T.D. ed, i.e., where the Alt est = 0.0 


Table 3.b) Planar/Non-conventional with 

Non-Complementary filter (r - 0.295 sec) 


* Nominal T.D. « 946' 


























V. 


Conclusions 


It was found that computationally configuration (lv) was simplest in 
concept. However, its severe sensitivity to the lateral y-deviation requires 
that this deviation be accounted for in precision usages of this configuration 
(Equation 17). A possible simplification is to use the average distance over 
the touchdown zone in conjunction with the azimuth measurement. 

On the other hand, configuration (1) proved to be the most cumbersome 
from a computational point of view. However, computing altitude from El^ and 
DME measurements only via the planar approximation of (Equation 11) , proved to 
be very accurate (within 'V 1/3 ft). Therefore, there is no difference compu- 
tationally if configuration (i) is used with the planar approximation and 
configuration (iv) with lateral y-deviation compensation using DME and Azimuth 
information. 

Configuration (iii) should yield similar or better results than con- 
figuration (1) if the same computation is used, since DME is used directly 
instead of the computed x distance. 

The nonconventlonal configurations require approximately twice the 
horizontal coverage as the conventional configurations. Furthermore, these 
configurations are constrained to place the antenna near the threshold. The 
nonconventlonal configurations require more than an order of magnitude greater 
vertical coverage compared to the conventional configurations. 

From the estimation point of view the combination of configuration (i) 
with the planar approximation and complementary filtering was best. (There- 
fore, configuration (iii) with the same filter should yield compatible or 
better results.) When configuration (lv) (Planar/nonconventional) was used 
with the y-deviation « 0 assumption and the complementary filtering, the 
results were not very favorable compared to the above case. However, when the 
lateral y-deviation is properly accounted for the estimation performance should 
be comparable with the above case. 

The noncomplementary (no acceleration aiding) filter for both cases was 
inferior in performance. The estimation errors were noisier due to the wider 
filter band-width compared to the complementary filter. 

Even though 20 ft (la) of DME noise did not show much affect in con- 
figuration (i) with either filter, we cannot give an acceptable DME noise 
value at this time. Neither can we give any conclusions on the impact of 
using the noncomplementary filter in a closed-loop (with the aircraft) flare 
auto-pilot. This point must be resolved with a more elaborate simulation of 
a controlled aircraft. 

The results are summarized in Table 4 below. 

Comparing configuration (i) to (iv) essentially results in a draw. 

There is little difference in computational requirements (for the (i) planar 
approximation and the (iv) exact equations) . Configuration (i) requires con- 
siderably less coverage chan (iv) ; however, better altitude and sink-rate 
estimates are available without complementary filtering with configuration (iv) . 



Table 4. Summary Table for El Analysis 
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